Growing static and dynamic length scales in a glass-forming liquid 
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We investigate the characteristic length scales associated with the glass transition phenomenon. 
By studying an atomic glass-forming liquid in negatively curved space, for which the local order is 
well identified and the amount of frustration opposing the spatial extension of this order is tunable, 
we provide insight into the structural origin of the main characteristics of the dynamics leading to 
glass formation. We find that the structural length and the correlation length characterizing the 
increasing heterogeneity of the dynamics grow together as temperature decreases. However, the 
system eventually enters a regime in which the former saturates as a result of frustration whereas 
dynamic correlations keep building up. 



The spectacular and universal character of the slow- 
down of relaxation in liquids as one approaches the glass 
transition seems to call for a detail-independent explana- 
tion that involves a growing length scale of one sort or 
another PQ[2]. Indeed, the issue of characteristic length 
scales has become one of the main themes of most the- 
ories of the glass transition. A large body of work has 
been triggered by the experimental and numerical finding 
that the dynamics becomes increasingly spatially hetero- 
geneous as it slows down [3]. Evidence for the existence 
of a related growing length scale has been obtained ei- 
ther from measures of various kinds of dynamic cluster- 
ing [SHS] or from multi-point space-time correlation func- 
tions [IHS]. On the other hand, despite the rather gen- 
eral proposition that a growing time scale should come 
with a growing static length scale [5] and evidence that 
the dynamical heterogeneities have some structural ori- 
gin [TU], the quest for relevant static correlations has 
been thwarted by obstacles. Among the latter is the well 
known fact that the simplest pair correlations probed by 
the static structure factor show no significant features as 
temperature decreases. 

There are several possible, and not mutually excluding, 
strategies to search for static correlations associated with 
the glassy properties of a supercooled liquid; such correla- 
tions presumably involve some "hidden" order parameter 
that does not show up in simple structural measures such 
as the static structure factor. One may then think of: 
(i) studying how far amorphous boundary conditions in- 
fluence the system, thereby looking at static point-to-set 
correlations pQ [HI [HI [H] , (ii) using finite-size analysis for 
various thermodynamic quantities [131 [H] , or else (iii) , 
provided that the locally preferred arrangement of the 
molecules in the liquid has been properly identified [T5l - 
[21], investigating the static pair correlations of the asso- 
ciated local order parameter, which amounts to consid- 
ering multi-particle correlations as in bond-orientational 
order parameters [H]. In this Letter, we follow the latter 
route. 

We investigate, at a microscopic level, the relationship 



between structure and dynamics in a glass-forming liquid. 
We identify static and dynamic lengths that both grow 
as relaxation becomes slower and the glass transition is 
approached. The model we focus on is motivated by the 
frustration-based approach of the glass transition [531 [53]: 
it is a 2-dimcnsional (2D) monodisperse 6-12 Lennard- 
Jones liquid whose specificity resides in the negatively 
curved (hyperbolic) space it is embedded in. This model 
in "flat" (Euclidean) 2D space is not frustrated because 
the local hexagonal order can freely propagate in space 
to form a triangular lattice. Placing the system in hy- 
perbolic geometry introduces frustration, whose strength 
is associated with the curvature of the hyperbolic plane. 
As a result the liquid does not crystallize and forms a 
glass upon cooling pi]. Frustrated hexagonal order in 
2D negatively curved space can be thought of as the 
analog of frustrated icosahedral order in 3D Euclidean 
space [22]. We have previously shown that the model 
displays the common features found in glassformers phe- 
nomenology [24] . 

In a nutshell, the hyperbolic plane is a 2D homoge- 
neous space of constant negative curvature —k^ with 
metric in polar coordinates expressed as d^s = d^r -|- 
(sinh(Kr)/K)^ d^Lp. The peculiarities of negatively curved 
space (boundary effects, parallel transport of vectors, ab- 
sence of global embedding in Euclidean space) entail sev- 
eral tricks for implementing molecular dynamics simula- 
tions. For details we refer to [24l - f26] . The simulations 
have been performed in the NVE ensemble at a density 
pa^ — 0.85 and for a frustration (curvature) parameter 
KCT from 0.02 to 0.2 (k = is the unfrustrated Euclidean 
plane). To enlarge the system size as needed to compute 
correlation lengths, we have considered periodic bound- 
ary conditions with both octagonal and 14-gonal elemen- 
tary cells. (In hyperbolic geometry, the elementary cell 
area A depends on the genus g of the associated com- 
pact manifold as, A — 4:TrK~^{g — 1), with 17 = 2 for the 
octagon and 3 for the 14-gone used here [25].) Over the 
range covered by the simulations, we are able to explore 
three different regimes in the liquid as temperature T de- 
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FIG. 1: Growth of hexatic local order: (a) Ge{r) 
normalized by its maximum value for na = 0.1 and 3 
temperatures: T/T* = 4.35, 1.19,0.59 (from bottom to 
top) . (b) Correlation or modulation length (extracted 
either from an exponential fit to the envelope of Ge{r) 
or from the location of the first modulation peak when 
present, see (a)) versus T* /T for 3 frustrations. 



creases |24j : the "normal" liquid range above a temper- 
ature T* ~ 0.75 that is characteristic of hexagonal-like 
ordering in the Euclidean space, the "fragile supercooled" 
liquid regime, which we have claimed to be controlled by 
the avoided transition at T* , and a regime dominated by 
the kinetics of rare irreducible frustration-induced topo- 
logical defects. The first two regimes presumably repre- 
sent the experimentally accessible range in fragile glass- 
forming liquids. 

Let us first discuss the static structure of the liquid. 
We have taken advantage of the fact that the locally 
preferred order is well identified in the present model: 
it is hexagonal, or more appropriately hexatic, for the 
small frustrations studied here. In addition to the stan- 
dard radial pair distribution function g{r), we have then 
monitored the extension of the local order through two 
kinds of observables: (i) the topological defects, whose 
local environment differs from 6-fold hexagonal arrange- 
ment (as obtained through a hyperbolic generalization 



of the Voronoi tessellation of particle systems), and (ii) 
the (hexatic) local bond orientational order parameter 
V'6(j) = (l/-^b) Xl<fc> ^^P(i6%fc) where the sum is over 
the N}, nearest neighbors of the jth atom and 9jk is the 
angle characterizing the "bond" between j and k. Its 
pair correlations are obtained from 
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where V'6(i|j) is the order parameter when parallel trans- 
ported from point to point rj along the path Tij, 
chosen here as a geodesic, and Vij is the distance be- 
tween the two points along Tij] the fraction with the 
delta function is appropriate for the hyperbolic metric. 
In addition, we have computed a measure of the fluc- 
tuations through X6 = -/V((*6*S) - (*6) (*5)), where 
= (1/iV) Y.^'<Pe{i\0) with the (arbitrary) origin. The 
results for GQ{r) at na — 0.1 are illustrated in Fig. la 
The envelope of the function decays at high T with a cor- 
relation length ^6 and as T is lowered, it shows a mod- 
ulation that is associated with the random pattern of 
frustration- limited hexatic domains (see Fig. 3 in [24j ) : 
as shown in the figure, the characteristic length then 
becomes a modulation length [37j . The T dependence of 
^6 is displayed in Fig. [lb] grows with decreasing T and 
saturates at a value ~ which, as seen from the con- 
figurations and their Voronoi tessellation [24], represents 
the average distance between the irreducible defects (neg- 
ative 7-fold disclinations) that destroy extended hexatic 
order. Note that due to the large system sizes required 
and the very time-consuming numerical procedures, we 
have not been able to compute the associated correlation 
length for the longest relaxation times and the lowest 
curvature (e.g., ^g at Ka = 0.05 has not yet saturated in 
our presently accessible range). 

The connection between spatial extension of the lo- 
cal order and slowdown of the relaxation is shown in 
Fig. [2] From the consideration of defect dynamics at 
low T [211 [2H] and a scaling hypothesis, we expect that 
the relaxation time r behaves as 
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and < 7/; < 1 [ig, and / (^e) ^ £.i exp{-E{^e)/T,) 
where is a low-T crossover below which hexatic or- 
der has saturated (see below). The best fit to the data 
gives an exponent ip ~ 0.5 (see Fig. [2]), but this value 
is only indicative due to the limited range of data. In 
connection with the previous observation that the super- 
Arrhenius character (fragility) increases with decreasing 
frustration [24j , this provides strong evidence that in the 
accessible T range, the slowing down of the dynamics is 
controlled by the (frustrated) extension of the local order 
due to the proximity of the avoided transition at T* |23j . 
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FIG. 2: Relation between the growth of the relaxation 

time T and that of the static length ■ scaled 
relaxation time £,q'^'^ log(Cff ^ t/too) (with (7=1) versus 



T* /T for the same frustrations as in Fig. lb The 
crosses correspond to the regime in which has 
saturated to . 



We now turn to the characterization of the spatial cor- 
relations in the dynamics. To this end, we characterize 
the average local dynamics and its fluctuations. From the 
distance travelled by any atom j during time i, dj(t), we 
compute the hyperbolic generalization of the local con- 
tribution to the "instantaneous" self intermediate scat- 
tering function, fsj{k,t) = P_ i^jfc (cosh(K(ij(i)), where 
P_i_i_jfc is a Legendre function of the first kind and k is 
chosen as usual near the maximum of the static structure 
factor [30]. From this, one has access to (i) the average 
dynamics Fs{k,t) = {^/N)J2j=i{fs,j{k,t)), from which 
we have extracted the relaxation time r used above 24J, 



(ii) the fluctuations XiW = {'^/N){iJ2Li Sfs,jik,t)f) 



with Sfs,j = fs,j - ifsj) 
tions in the dynamics, 



and (iii) the spatial correla- 
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whose integral over the whole system area gives back 
X4{t)- (As for Eq. (1), Eq. (2) reduces to the familiar 
Euclidean expression when k — > 0.) We have focused on 
the maximum of X4(^) ^-^^d on G4{r,t) for the time 

tmax ~ T at which this maximum occurs; from an expo- 
nential fit of the envelope of the latter function, we have 
extracted a dynamic correlation length ^4. 

A comparison between the "dynamic" length scale ^4 
and the static one is displayed in Fig. [3] The T depen- 
dences of the two quantities are strikingly similar, with 
both lengths saturating at low enough T (when accessi- 
ble) to a value determined by the curvature (^g roughly 
saturates to and ^4, which is shifted and divided by 
a T-independent factor in Fig. [3j to about 0.5k~^); ^4 
and ^g appear linearly related to each other. This shows 
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FIG. 3: Static (,^g) and dynamic (^4) length scales 
versus T* /T for 3 frustrations: from bottom to top, 
Ka = 0.2,0.1,0.05. ^4 is shifted and divided by a 
T-independent factor. 



that the buildup of spatial correlations in the dynam- 
ics, which is associated with growing dynamical hetero- 
geneities, and the one in the static local order parame- 
ter, which has been seen above to control the slowdown 
of relaxation, are directly connected in the T range from 
normal liquid to the beginning of the regime dominated 
by irreducible defects. 

Although we only have limited information on the low 
T behavior below some crossover T^,, once the lengths ^g 
and ^4 have saturated, it is interesting to push further the 
analysis of this regime. We anticipate that hexatic order 
no longer evolves because it has reached its frustration- 
induced limit. This is confirmed by the behavior of the 
static "susceptibility" XeiT) (not shown here) that sat- 
urates just as ^g does. On the other hand, the dynamic 
susceptibility xT"'^ keeps growing with decreasing T, al- 
beit with a change of rate: xT""^ roughly crosses over to 
a 1 /T'^ behavior (see Fig. 4) . Spatial correlations are still 
increasing as the dynamics slows down in the irreducible- 
defect regime. A close examination of t/4(r, tmax) hints at 
what may happen: when ^4 saturates, another length ^4 
emerges, which in our model corresponds to a faster than 
exponential decay of the envelope of Q4 beyond the sat- 



uration length (see Fig. 4a I. Although we cannot follow 
the growth of this new dynamic length, we can indirectly 
measure its effect by comparing, as illustrated by Fig. |4b[ 
the dynamic susceptibility xT""^ obtained for the full sys- 
tem area and that computed by truncating the spatial in- 
tegration of G4:ir,tmax) to a cutoff distance Tc ~ 0.8k~^, 
i.e. larger than ^4 at saturation but smaller than ^4. 
The latter saturates as T decreases whereas, as stated 
before, the former keeps increasing with just a change of 
rate. We can thus conclude that in the irreducible-defect 
regime below T^,, there is a decoupling between the spa- 
tial correlations in the dynamics and in the frustrated 
hexatic order. To our knowledge this is the first time 
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FIG. 4: Crossover in the space-time correlations: (a) 
Giir^t^ax) for T/T* = 0.82 and na = 0.2; ^4, 
characteristic of the exponential decay, is near 
saturation but an additional length ^4 > ^4 emerges, 
(b) For the same conditions, comparison between xT""^ 
obtained for the whole system area and by integrating 

QA{r,tmax) up tO Tc ^ 0.8k"^ 



that such a phenomenon has been observed. 

We stress that the crossover to the irreducible-defect 
regime is likely to be intrinsic to frustration rather than 
a mere peculiarity of hyperbolic space. In any system in 
which growth of the locally preferred order is frustrated, 
one expects that there is an intrinsic frustration length at 
which the spatial extent of the local order saturates |23J . 
For instance, this should be the case in the 2D models 
of spherical polydisperse particles when polydispersity is 
higher than the threshold value above which crystalliza- 
tion in hexagonal structures is avoided [17] , as well as in 
3D models with icosahedral local order [131 UHl HH] • In 
the temperature range in which saturation has occurred, 
there is an irreducible density of defects that essentially 
stays constant with decreasing temperature. If frustra- 
tion is not too strong, so that the average distance be- 
tween defects is large, the initial part of this irreducible- 
defect regime is expected to lead to simple Arrhenius 
dependence of the relaxation time, with the activation 



barrier determined by the energy scale of the defects. 
What happens at still lower T is only subject to specu- 
lation. In principle, the defects at least slightly interact 
and may be subject to kinetic constraints induced by the 
locally ordered environment. 

The present study provides insight into the structural 
origin of the main characteristics of the dynamics (slow- 
down of relaxation and increasing heterogeneous charac- 
ter) in a glass-forming liquid. Interestingly, it is found 
that frustration, a phenomenon which has been hypoth- 
esized to be ubiquitous in liquids [TTl [22l [23] , imposes a 
saturation of the spatial extent of the local order lead- 
ing to a crossover to a low-T regime dominated by rare 
irreducible topological defects. The slowdown and het- 
erogeneous character of the dynamics then decouple from 
the frustrated local order. Whether this further results in 
a simple noncooperative regime or involves a new type of 
static correlations remains to be investigated. A first step 
in this direction would be to compute defect-defect corre- 
lations as well as point-to-set correlations in the present 
model. 
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